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ABSTRACT 

Aims. In some of the lensed quasars, color differences between multiple images are observed at optical/near-infrared wavelengths. 
There are three possible origins of the color differences: intrinsic variabilities of quasars, differential dust extinction, and quasar 
microlensing. We examine how these three possible scenarios can reproduce the observed chromaticity. 

Methods. We evaluate how much color difference between multiple images can be reproduced by the above three possible scenarios 
with realistic models; (i) an empirical relation for intrinsic variabilities of quasars, (ii) empirical relations for dust extinction and 
theoretically predicted inhomogeneity in galaxies, or (iii) a theoretical model for quasar accretion disks and magnification patterns in 
the vicinity of caustics. 

Results. We find that intrinsic variabilities of quasars cannot be a dominant source responsible for observed chromatic features in 
multiple quasars. In contrast, either dust extinction or quasar microlensing can nicely reproduce the observed color differences between 
multiple images in most of the lensed quasars. Taking into account the time interval between observations at different wavebands in 
our estimations, quasar microlensing is a more realistic scenario to reproduce the observed color differences than dust extinction. 
All the observed color differences presented in this paper can be explained by a combination of these two effects, but monitoring 
observations at multiple wavebands are necessary to disentangle these. 

Key words. Accretion, accretion disks - h Gravitational tensing - dust, extinction - quasars: general 

1 . Introduction jects is important to obtain deeper insight into the effect of chro- 

maticity in lensed objects. 



Gravitational lensing represents a major observational tool in An important example appears as the color difference Q be- 

H modern astronomy and astrophysics. Although its practical ap- tween multiple images of lensed quasars. Although multiple im- 

^ plication has been mainly focused on cosmological aspects at ages of lensed quasars should have the same color as the non- 

the beginning (e.g.,|Zer dovjcE 1 1 964 |Refsda| 1 1 964 [Paczynskil, lensed image according to the principle of gravitational lensing, 

1986), the subjects extend to studies of galactic structure (e.g., some lensed quasars show a clear color difference between dif- 

Oguri et alj 12002), searches for exo-planets including earth fgrent images. One possible explanation for the observed chro- 



mass planets (e.g., [BeauUeu et aj], |2o66i) and others. A fun- maticity is differential dust extinction inside the lens galaxy. It 
damental feature of gravitational lensing is achromaticity; this is apparent that the spatial distribution of interstellar gas and 
means that in principle gravitational lensing effects show no de- dust is inhomogeneous in galaxies. The Ught path of different 
pendence on wavelength. In many situations, this achromatic- images pass different parts of the lens galaxy, and it is natu- 
ity is important to discriminate gravitational lensing effects from ral to consider that different images are affected by spatially 
other annoying phenomena for the detection; for instance, vari- varying dust extinction characteristics. Within this scenario, the 
able stars in galactic microlensing surveys (e.g., |AIcock et al.t amount of color difference between images at given wavebands 
11996 '). However, some gravitational lensing phenomena are as- is mainly determined by three quantities; the redshift of the lens 
sociated with unexpected chromatic feature s; that is, observe d galaxy, the column density of dust, and the extinction proper- 
properties have wavelength dependence (e.g.,|FaicoietaU[i99^ ties of dust. Thus, we can probe these three quai itities from 
In such cases, the gravitational lensing hypothesis has been con- photometric data of more than three wavebands. iFaIco et al.l 
firmed using other supporting observations; studying these ob- (1999) have investigated a method to probe the dust extinc- 
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' Throughout this paper, we use the word 'color difference' as differ- 
ence between colors of multiple images, and never use it as difference 
between different colors. 
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tion in distant g a laxies , iToft et alj f2000l) : iMunoz et alj ( |2004 : 
lEliasdottir et alj (l2006h : iMediavilla et all ( 12005 ') have explored 
the extinction law at high-redshift galaxies, and Dai et al. (200^ 
have probed the dust-to-gas ratio of galaxies at cosmological dis- 
tance by combining the hydrogen column density obtained with 
th e Chandra X-ray Observatory and the color excess obtained 
by 'Falco et alJ (1199 9). This enables us to measure the absorption 
properties at distant galaxies directly, and it is complement to 
emission- weig hted measurements of dust properties. In contrast, 
iJean & Surdeil (11998) have proposed a method to estimate the 
redshift of the lens galaxy under an assumption on the extinction 
properties. Direct detection of the emission from the lens galaxy 
is not necessary for this method, and it can be a strong tool to 
estimate redshift of faint or so-called "dark" lens galaxies. 

Unfortunately, differential extinction is not the only scenario 
able to explain the observed chromaticity in multiple quasars. 
Because multiple images have different light paths, a delay in the 
arrival time is always present between the images. This indicates 
that we observe slightly different epochs of the lensed quasar at 
the same time via multiple images. Quasars intrinsically change 
not only their luminosity but also their color with time (e.g., 
IWilhite et alj,l2005l) . Incorporating time delay and intrinsic color 
variabilities, one could explain the observed chromaticity to 
some degree. In addition, quasar microlensing can be yet another 
candidate to explain the observed chromaticity. A theoretical 
investigation of chromatic f eatures of quasar microlens ing has 
initially been made by Wa mbsganss & Paczviiskil (Il99ll) in the 
case of so-called "Huchra's lens (Q2237-H0305)". In other sys- 
tems, possible evidence for ch romaticity due to mic r olensing has 
been r eported, for example, bv Burud et al.l ( l2002al) : rNakos et al] 
(I2005 h. Since there is no essential difference between Huchra's 
lens and other multiple quasars, color changes due to quasar mi- 
crolensing can be an origin of the observed chromaticity and can 
work as a contaminant in exploring dust extinction in the lens 
galaxy. 

In summary, there are three possibilities for the chromaticity 
in a multiple quasar: (i) intrinsic quasar variability, (ii) differ- 
ential dust extinction, and (iii) quasar microlensing. In this pa- 
per, we consider all these three possible scenarios to explain the 
observed chromaticity with realistic theoretical models and reli- 
able empirical relations. We will briefly introduce the observed 
chromaticity in section 2. The three possible scenarios for the 
observed chromaticity are individually examined in sections 3, 
4, and 5. Discussions about possible and realistic origin of the 
observed chromaticity are presented in section 6. Throughout 
this paper, we adopt the following cosmologic al parameters : 
7 0.3, Qa = 0.7, and Hq = 70 kms 'Mpc-' dSpergel et all 
120061) . 



2. Observed Chromaticity 

Several observational programs for lensed quasars are currently 
carried out, and we here use the data taken by one of such pro- 
grams, CASTLES 0. CASTLES provides us with photometric 
data of lensed quasars of high and equal quality, and the data are 
suitable for comparison with theoretical models and predictions. 
Roughly 100 lensed quasars have been listed so far, but we pick 
up only 25 objects whose photometric data at F160W of Hubble 
Space Telescope (HST) NICMOS, and at F555W and F814W of 
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^ CfA-Arizona Space Telescope LEns Survey. The 
data and details of this program are presented in 
|http : // www . cf a . harv ard . edu/castle s/ See also publica- 
tions of CASLTES fe.g.. lLehareTaill2000l) for the details. 
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lens redshift (zj 

Fig. 1. The lens redshift (z/) and the source redshift (z,) of the 
current samples are presented by the filled circles in this fig- 
ure. Objects without the lens redshift are indicated by the fiUed 
squares with arrows. 



WFPC2 are publically available online 0. This sample en- 
ables us to use two independent colors to disentangle the origin 
of the observed chromaticity. This is useful for breaking some 
degeneracy which will easily occur from a single color informa- 
tion, for example, degeneracy between extinction properties of 
dust and amount of dust. Details of the sample lensed quasars 
and lens galaxies are summarized in Tables[Tl and|2] 

Here, we briefly summ arize how t h e pho tometric data are 
analyzed and obtained by Lehar et aTl ([2000). The brightness 
of objects was measured by minimization among the ob- 
served images and models for lensed quasars and lens galaxies. 
A lensed quasar image is modeled by a PSF with a sharp central 
core (FWHM ~ 0.1 arcsec) and an outer skirt of low intensity 
(< 1% of the peak intensity), and lens galaxies are modeled by 
an ellipsoidal exponential disk or de Vaucouleurs models. Each 
lensed quasar image has 3 parameters for fitting (2-dimensional 
position and flux), and each lens has 6 parameters for fitting (2- 
dimensional position, flux, major axis, axis ratio, position angle 
of major axis). Uncertainties of these decomposition procedures 
are included in errors presented in Tables [1] and|2] 

The redshifts of the lens (z/) and the source (z,) listed 
in Tables [1] are also taken from CASTLES, and presented in 
Figure [1] The redshifts of all the lensed quasars are known as 
presented in Tables [T] In contrast, the redshifts of the lens galax- 
ies are not known in all cases. Objects with unknown redshifts 
or with tentative redshifts are listed at the bottom of Table [1] For 
these objects, the upper limits of their lens redshifts are set to be 
the source redshifts and are denoted by the filled squares with 
arrows in Figure [1] The source quasars and the lens galaxies are 
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Table 2. Summary of lens galaxies in the current sample. Their magnitude and error at F160W, at F555W and at F814W are pre- 
sented. All the data are taken from CASTLES web page. Types of the lens galaxies indicated by their spectrum, color or brightness 
profile are also shown. Early and late type is presented by 'E' and 'L', respectively. For a sample with two lens galaxies, they a re de- 
noted by Gand G'. Referenc es. (l)|Lehar et al. 2000; (2) Tonrv & Kochanek 1999; (3) Fassnacht & Cohen 1998; (4) Keeto n et al. 



19981; (5) iBurud etaLl]l998l; (6) iLubin et al. 2000; (7) Younsetal. 1980; (8) Angonin-Willaime e t al..,1994; (9) Kochane k et al 



200C ; ( 1 0) iRemv et'af 1 1 998t ( 1 1 ) Impey et al. 1998; (12) Impey et al. 1996; (13) Chavushyan et al.l ll997l; (14)|Koopmans & Treu 



2002L r 1 5) lFassnachtetanfT999t (16 VlLopez etal.lll998t (17) lHuchra etal.lll985L 
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distributed at z ^ \-4 and at z ^ 0-1, respectively 0. In total, now 
we have 25 objects (or 40 image pairs) for our current purpose. 

In TablefT] the colors of each lensed quasar image for F555 W 
- F160W and F814W - F160W, and the F160W magnitudes of 
each image are also presented together with their errors. The er- 
rors are evaluated in the usual manner, namely by calculating the 
square root of the sum of the squares of the individual error. It 
is apparent from Table [1] that some of the lensed quasars show 
color differences between different images. This feature can also 
be found more clearly in Figure |2j the color differences of mul- 
tiple images relative to that of the brightest image at F160W, 
A(wf555iv -ffifieoiy) and A(mf8i4H' -'«fi60w)- Images located at 
the lower left part in Figure|2]are relatively bluer than the bright- 
est image at F160W in the system. Gravitational lensing effect 
has no wavelength dependence in principle, and the magnifica- 
tion factor should be the same at different wavebands. Since the 
magnitude differences between images correspond to the ratio 
of magnification factors only if the macro lensing contributes 
to the magnification, the magnitude difference should be iden- 
tical among any wavebands. That is, observational data points 
should be located at the origin, (0, 0), in Figure |2] within the er- 



lEigenbrod e t al. (2007) has recently measured that the lens red- 
shifts of FBQ0951+2635 and HE2149-2745 are 0.260 and 0.603, re- 
spectively. However, these redshifts are within this redshift range of the 
lens galaxies and our conclusion are not modified by the measurements. 



ror. Though some of the multiple quasars show achromaticity as 
expected, others show anomalous chromaticity; i.e., some data 
points in Figure |2] significantly deviate from the origin, ranging 
up to ~ 2 mag. 

Although there is a wavelength overlap among these two 
combinations of the color differences, F555W - F160W and 
F814W - F160W, we choose these combinations because of 
their small errors compared to combinations of F555W - 
F8 14W. The lower-left to upper-right trend of observational data 
points in Figure |2] or following similar figures is a natural con- 
sequence of the wavelength overlap, and hereafter we just focus 
on deviations of observational data from the expected values of 
our calculations rather than the trend. 

Here, we review what kind of effect should be included for 
interpretation of the observed chromaticity. By using the intrin- 
sic magnitude of a quasar (mint), the expected magnitude of an 
image of multiple quasars (mobs) at wavelength (/l) at any given 
time (f) is written as follows: 

mobs(^, i) = mint(/l, f - t) + A |y7^) ~ 0] , (1) 

where t. A, and yu indicate the arrival time delay due to gravi- 
tational lensing, the extinction along a path from the source to 
observer, and the magnification due to gravitational lensing, re- 
spectively. To include the finite source size effect, as is differ- 



4 



Atsunori Yonehara et al.: Origin of Chromatic Features in Multiple Quasars 




A(mp 



2 

J [mag; 



<l 



(b) 




-2 - 



A(mj 







[mag] 



Fig. 2. Color differences between multiple images are plotted with error bars. For each system, the brightest image in F160W has 
been used as color reference. The abscissa and ordinate represent the color difference derived from F555W and F160W {t!^{mfsisw - 
mpi^Qw)), and that from F8 14W and F160W (A(OT/rsi4H' -'Wfieow)), respectively. Panel (a) is for lensed quasars for which the redshift 
of the lens galaxy is known, while Panel (b) is for lensed quasars for which the redshift of the lens galaxy is unknown. 



ent from a usual expression, the magnification due to gravita- 
tional lensing is presented as a function of wavelength and time 
in equation[T](see section 5 for more details). These three terms 
all affect the multiple images with different magnitudes for dif- 
ferent images. Thus, they can be possible origins of the observed 
chromaticity. In the following three sections, we examine how 
these effects may produce the observed chromaticity. 

3. Intrinsic Variabilities of Quasars 

It is well known that quasars show temporal flux variabilities at 
any waveband s, and the color becomes bluer when the flu x gets 
brighter (e.g., IWilhite et all 120051: ICristiani et all 11997 ). This 
means that the color of quasars changes with time according to 
the flux variabilities. Additionally, the arrival time delay due to 
gravitational lensing is different at different images. This pro- 
duces the relative time delay between multiple images of the 
same quasar, and multiple images observed at the same epoch 
correspond to images at intrinsically different epochs. Therefore, 
it is probable that multiple images observed at the same time 
show different colors. To estimate this effect, we evaluate typical 
time delays between images and typical color change of quasars 
for a given time interval. 

Though the actual time delay is determined by the density 
profile of the lens galaxy, the location of the source respective to 
the lens on the sky, the source redshift, and the lens redshift, we 
can roughly evaluate the typical time delay from observed prop- 
erties of multiple quasars. By assuming that the density profile 
of a lens g alaxy is approximated b y a singular isothermal sphere 
(SIS; e.g.. ISchneider et akl Il992h . the time delay between two 
images, one image with positive parity (located at 9+ from the 
lens galaxy) and another image with negative parity (located at 
0- from the lens galaxy), is expressed by 

_ 1 + Zl DqiDqs 2\ 



.30(ilii)(^^^)(^)days, (2) 
\ 2 j\ 1 Gpc / \ 1 ai-csec2 / ^ 

where D,,!, Dos, and Dis represent the angular diameter distances 
from observer to the lens, from observer to the source, and from 
the lens to the source, respectively. Since the image separation is 
on the order of 1 arcsec in most multiple quasars, we can roughly 
estimate 0^ - 6*^ to be ~ 1 arcsec^. Thus, we adopt 30 days for a 
typical time delay between multiple images. 

Further, by using recent studies of intrinsic variabilities of 
quasars, we can also evaluate the color change of quasars on any 
given timescale. IVanden Berk et alj ( l2004l) : llvezid et all (120041) 
have investigated how quasar variabiUties in the rest-frame opti- 
cal/UV regime depend on other observational quantities such as 
rest-frame time lag. They analyzed relations among these quan- 
tities and the structure funct ion (a commonly used statistical 
measure for variabilities; e.g.. lKawaguchi et allll998 ) by using 
imaging data of quasars obtained by Sloan Digital Sky Survey. 
The resulting structure function (V) is expressed as 

V^(l+ 0.024M,) — mag, (3) 

\ ^RF j 

where M,-, AfRp, and /Irf represent the / band absolute magni- 
tude of a quasar in units of magnitude, the rest-frame time lag 
between observations in units of days, and the rest-frame wave- 
length in units of A, respectively. Equation [3] indicates that a 
larger flux change of quasars occurs at fainter quasars, at longer 
time lag, and/or at shorter wavelength. The absolute / band mag- 

nitude of the q uasar sampl e is distributed from ~ -21 to 30 

dVanden Berk et all |2004|) . Since Af^p and /Irp have the same 
redshift dependence, both of these quantities can be replaced by 
the observer-frame quantities. 

Substituting the typical time delay above (30 days) and 
the effective wavelength of each filter into equation [3] we ob- 
tain the expected flux variation of quasars. The expected val- 
ues for magnitude difference {6m) with a time lag of 30 days is 
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0.043 - 0.075 mag, 0.059 - 0.105 mag, and 0.053 - 0.093 mag 
at F160W (mean wavelength: 16071 A), F555W (mean wave- 
length: 5337 A), and F814W (mean wavelength: 7900 A), re- 
spectively. Here, the maximum and the minimum value of mag- 
nitude differences are obtained for M, - -21 and M, - -30, 
respectively. Even if we take into account filter responses, these 
values would not change much. The wavelength coverage of all 
the filters is less than ~ +20% of the wavelength center (see also 
the upper panel of Figure|4|, and the expected flux variations are 
less than +7% from the above values within the coverage (see 
equation[3]|. By comparing these expected values for magnitude 
difference at different wavebands, we are able to estimate the 
possible range of the expected color differences. 

In general, the range of the expected color difference be- 
tween two images for two bands labeled with j and k, \A(mj - 
mk)\, is estimated by \\6mj\ - \6mit\\ < \A(mj - mk)\ < \6mj\ + \Smit\. 
The maximum value corresponds to the case where variabilities 
in band j and that in band k have completely negative corre- 
lation, i.e., when a quasar gets brighter in band j, a quasar al- 
ways gets fainter in a A; band. The minimum value corresponds 
to the case where variabilities in these two bands have com- 
pletely positive correlation, i.e., when a quasar becomes brighter 
in band /', it always gets brighter also in band k. For quas ars with 
-21, the fainter end of the sample of Vanden Ber k et al] 
12004)), the expected color differences are ranging from |0.105 - 
0.0751 = 0.030 mag to 0.105 + 0.075 = 0.180 mag for 
|A(mF555w - »Jfi60w)I, and from |0.093 - 0.075| = 0.018 mag 
to 0.093 + 0.075 = 0.168 mag for |A(mF8i4w - '«fi60w)I- Since 
the observed chromaticities as shown in Figure |2]have values up 
to ~ 2 mag, it is impossible to reproduce all the observed chro- 
maticity only by this scenario which is able to reproduce a color 
difference up to 0.2 mag. For intrinsically brighter quasars, the 
expected flux change within a time interval is rather small, and 
the expected color differences are smaller than the observed val- 
ues. Moreover, quasar variabilities at any two wavebands pre- 
sumably have positive correlat ion with each o ther rather than 
a negative correlation (e.g., Wil hite et al.L l2005h . and the actual 
color difference should be smaller than our estimated maximum 
values. 

Of course, such positive correlation will be lost if ob- 
servations for an object at different wavebands are carried 
out at different times. Denoting this observational interval 
among two bands (labeled with j and k) as f/oj, the color 
difference including time delay (r) is expressed as A(mj - 

m) = [mj{t)-mkit + ti„g)j - \mj(t + T)-mk{t + T + ti„g)j = 

nijit) - nijit + T)j-|^m^(f + tiag) - nikit + t + tiag)^. This relation 
indicates that the expected color differences are a simple com- 
bination of magnitude differences at each band adopted above. 
Since an actual time delay is nothing to do with the observa- 
tional interval among different wavebands, tiag is the only factor 
that makes difference from the previous sit uation, i.e., withou t 
the observational interval. As presented bv iLehar et alj (l2000l) . 
the observational interval among different wavebands spanned 
up to tiag ~ 2 years for some objects, and it is comparable 
to m aximum timescale of quasar variabilities (e.g..' lvezic et al.L 
l2004) Fl Since physical correlation could be lost after a timescale 



^ Generally, structure function of intrinsic quasar variability con- 
sist of two parts; power-law component below a certain timescale 
and flat component above the timescale. Since the flat component 
of the structure function originates from only a random process, the 
timescale where the power-law component ends corresponds to maxi- 
mum timescale of quasar variability due to some physical process. 



longer than maximum timescale of quasar variability, a com- 
pletely negative correlation among two wavebands is even pos- 
sible when the observational interval is comparable to or longer 
than ~ 2 years. However, even if it is the case, the amplitude 
of variations estimated above fixes upper bounds on the ex- 
pected color differences, and it is impossible to produce color 
differences larger than |A(mF555w - 'wfi60w)I - 0.180 mag and 
|A(mF8i4w - '«Fi60w)l = 0.168 mag. Thus, the expected color 
differences should be still within the range estimated above, and 
thus the observed color differences cannot be explained by in- 
trinsic variabilities of quasars alone. 

The model for the lens galaxies that we used here for esti- 
mating the time delay between multiple images is rather sim- 
ple, though the time delay depends on the applied lens model. 
This has a lready been mentioned in previous studies (e.g., 
lO guri et aLl [2002), and the expected time delay of more real- 
istic lens models can be reduced down to an order of magnitude 
from that of SIS. Diversity of the expect ed time delay is al so in- 
vestigated by different approaches (e.g., ISaha et all l2006l) . and 
the time delay range would be typically between ~ 10 and ~ 300 
days. Further, if the lens redshift is unknown, the expected time 
delay of any lens model can have a dispersion of one order of 
magnitude for a given image separation. Consequently, taking a 
more realistic lens model into account, the expected time delay 
can change up to two orders of magnitude. Even if this is the 
case, the expected color difference between multiple images can 
change only by a factor of 4 (see equation|3]l. 

From the observational point of view, the measured time 
delay shows wide variety as predicted by the previous theo- 
retical studies (e.g., lOguri et al.L 120021; ISaha et al.L 12006 ). and 
the measured values are different from 30 days. In Table [3] 
we summarize the time delay between multiple images in 
the current sample. The time delay of objects which are not 
listed in Table [3] Q0142-100, MG0414-H0534, B0712+472, 
SBS0909+523, LBQS 1009-0252, B1030+071, MG2016+112, 
B2045+265, Q2237+0305, QJ0158-4325, APM08279-H5255, 
BRI0952-0115, Q1017-207, Q1208+101, and FBQ1633H-3134, 
has not been successfully measured yet. As we can see in 
Table [3] some multiple quasars have one order of magnitude 
lon ger time delay than the typical value applied here (e.g., 
Kundic et al.L 1 19971) . In such systems, of course, the expected 
flux change is larger than that we estimated for multiple quasars 
with the time delay of 30 days. However, the color difference 
becomes only twice the value estimated above (see equation[3]l, 
and the observed color differences are not still explained. Again, 
it is clear that intrinsic variabilities of quasars alone cannot re- 
produce observed color difference between multiple images. 

4. Differential Dust Extinction 

The light paths of different images intersect different positions 
of the lens galaxy. In general, the column density of gas and dust 
along different lines of sight differs in all galaxies, and differ- 
ent images experience different levels of dust extinction. We can 
also consider that the shape of extinction curve is different at dif- 
ferent positions in the lens galaxy. For simplicity, however, we 
basically take into account only the inhomogeneity of dust col- 
umn densities, but we also examine the difference among various 
empirical extinction curves derived from the Milky Way (MW) 
and the Smafl Mageflanic Cloud (SMC). 

We here apply an extinction law derived for the MW dust 
by Cardelli et al. ( 1989). In this extinction law, dust extinction is 
characterized basically by two parameters; Ay, dust extinction at 
y band which reflects amount of dust, andRy [= Av/Eb-v (Eb-v 
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Table 3. Time delay between multiple images in current sample is presented in units of day. Wavebands used for t he measure- 
ment are also presented with the reference. Only objects with the available time delay are listed here. References. (1) Biggs et alj 
19991: (2) Hiorthetal. 2002; (3) UU an et al. 2006; (4) Kundic et al. 1997; (5) Haarsma et al. 1999; (6) Ofek & Maoz 20 03; 



Schechteret al.ll 19 97: (8) P atnaik & Narasimhall2001l: r9') lBurud et alJl2002bt riO') iKoopmans et al.ii200dl: Vin lBurud et aLli2002; 



(12) lJakobsson et al..i2005> 



Object and Image Pair 


Optical Delay 
(day) 


Waveband 


Radio Delay 
(day) 


Waveband 


B0218+357 A-B 






10.5 ±0.4 


8.4 and 15 GHz (1) 


RXJ0911+0551 A-B 


146 ±8 


I-band (2) 


— 




SBS0909+523 A-B 


45+' 


R-band (3) 






Q0957+561 A-B 


417+^ 
420+^ 


g-band (4) 
r-band (4) 


397+15 


6 cm (5) 
4 cm (5) 


HE1104-1805 A-B 


161+? 


V- and R-band (6) 






PGl 115+080 A-C 


9.4 


V-band (7) 






C-B 


23.7 ± 3.4 


V-band (7) 






B 1422+231 A-B 






1.5 + 1.4 


8.4 and 15 GHz (8) 


A-C 






7.6 + 2.5 


8.4 and 15 GHz (8) 


B-C 






8.2 ± 2.0 


8.4 and 15 GHz (8) 


SBS 1520+530 A-B 


130 ±3 


R-band (9) 






B 1600+434 A-B 






47!^ 


8.5 GHz (10) 


HE2 149-2745 A-B 


103 ± 12 


V- and i-band (11) 






FBQ095 1+2635 A-B 


16±2 


R-band (12) 





is the B -V color excess)]. Since Ay is physically related to the 
column density of dust, we can estimate the differential dust ex- 
tinction from differential column densities of gas by assuming 
values of the dust-to-gas ratio and Ry- In order to obtain differen- 
tial column densities of gas, we utilize a 2-dimen sional galactic- 
scale hydrodynamical simulation presented by 'Hirashita et al.l 



(I2003h . who adopt the calculation code of Wada & Norman 
(1200 Ih . The result clearly shows a clumpy distribution of gas 
in the simulated galaxy. The cumulative probability distribution 
of differential column densities, AA^h, is shown in Figure [3j we 
randomly select lO*" pairs of columns on the 2-dimensional grid 
points in the simulated galaxy, and derive the probability dis- 
tribution function of the difference of gas column densities be- 
tween each pair0. As we can see from Figure[3] the expected dif- 
ferential column density from the simulation spans a large range: 
The median value is ~ 3 x 10^*^ cm"^, and 90% of data is dis- 
tributed from ~ 6 X 10'^ cm"^ to ~ 7 x 10^' cm" ^. We apply the 
MW dust-to-gas ratio provided bv Bohlin et al.l ( 1197 8), and ob- 
tain Av ^Rvx Eb-v = 3.1 X A^h/(5.8 x lO^') = 5.3 x IQ-^^Nn, 
where A^h is the hydrogen column density in cm"^. In this pa- 
per, we apply the above relation also when we adopt values for 
Ry other than 3.1. We then can convert the expected differential 
column densities into the expected differential extinction in the 
V band, AAy. This conversion is fairly simple and the scale for 
this value is also denoted in Figure[3] The corresponding values 
of AAy for the median, 5%, and 95% of Figure [3] are ~ 0.1 mag, 
~ 3 X lO"-' mag, and ~ 3 mag, respectively. 

Additionally, it is well known that the SMC exhibits a very 
different extinction curve compared to the MW curve. The SMC 
curve lacks the 2175 A bump and rises steeply toward shorter 
wavelengths in the UV. We also adopt an SMC extinction curve 
in the current study. Several types of the SMC extinction curves 
have been provided, but we investigate the extinc tion law by us- 
ing an observational data of iGordon et al.l ( l2003 h (see appendix 
for details) which covers a relatively wide wavelength range. 



* A lensed image is usually located at the opposite direction from 
the other image with respect to the center of the lens galaxy. Since the 
gravitational evolution of clumps in galaxies is coherent only within a 
Jeans scale, such coherency may have a negligible effect on our result, 
i.e., spatial inhomogeneity over the galactic scale. 



AAy [mag] 
10-^ 10-2 IQ-l 



10^ 




10 



IQ19 10=^° 10^1 



10^ 



102 



Fig. 3. Cumulative probability distribution of differen tial column 
densities of gas obtained from lHirashita et a n (l2003h . The lower 
abscissa indicates the differential column density of gas in units 
of cm"^. The corresponding differential extinction at V band (in 
mag) is also denoted at the upper abscissa. Cumulative probabil- 
ities of 5%, 50% and 95% are indicated by the dashed lines. 



Some examples of the extinction cu rves that we adopted i n this 
study are shown in Figure |4] Since ICardelli et al.l (Il989l) have 
presented the extinction law only for Ry = 2.6 - 5.6, and since 
there is no guarantee that the empirical extinction law is valid 
for j^v < 2.6or j^v > 5.6, we put in this paper the same limit as 
ICardelUet all([T989l) . 

Denoting the response of the filter k and the spectrum of the 
source as fk{A) and S(A), respectively, the expected magnitude 
at filter k, nik, including dust extinction is expressed as 



nikiNn) = mi,„,k-2.5 log 



10 



-0.4A[A'h„V(1+::;)] 



MA)S(A)dA 



f^'-MA)S(A)dA 



.(4) 
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15000 
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A [A] 

Fig. 4. Extinction curves for several parameter sets (lower panel) 
and filter responses (upper panel) are presented. In the lower 
panel, the solid, the dotted, the dashed, and the long-dashed 
lines correspond to the extinction curve of the SMC with A^h - 
10^' cm-2, that of the MW with A^h = 10^' cm^^ and Rv = 5.6, 
that of the MW with A^h = 10^' cm"^ and /Jv = 3.1, and that of 
the MW with A^h = 3 x 10^' cm"^ and Ry - 3.1, respectively. In 
the upper panel, filter responses for F555W, F8 14W, and F160W 
are presented by the dotted, the dashed, and the solid lines, re- 
spectively. The abscissa of the lower panel is the wavelength at 
the rest frame of the absorber (or the lens galaxy) in units of A, 
and that of upper panel is the wavelength at observer's frame. 



we can estimate the expected color differences due to differen- 
tial dust extinction. The results are presented in Figure |5] 

As clearly seen in Figure |5] the differential dust extinction 
of 95% level (see Figure |3j can produce up to ~ 10 mag of 
color difference between multiple images of lensed quasars 0. 
The expected color difference is large enough to reproduce the 
observed chromaticity. Moreover, the slope nicely reproduces 
the observational data. We note that the above theoretical pre- 
dictions do not depend strongly on the lens redshift and the ap- 
plied extinction laws (see Figure|5j. However, the expected color 
difference can be as large as a few magnitudes, exceeding the 
observed values for the color difference. This may be because 
we have adopted the simulation of a relatively gas-rich galaxy 
that has n ot yet converted a substantial fraction of the gas mass 
into stars dHirashita et al.L l2003t) . Actually, many lens galaxies 
are suggested to be early type galaxies, and the dust amount and 
distribution in such galaxies must differ from what we have ap- 
plied here. In this respect, our approach may be too simplistic. 

In the lower panels of Figure |5] some of the expected color 
differences are presented together with the same observational 
data as those in Figure |2] for comparison. Most of the observa- 
tional data in both of these panels (panels c and d are respec- 
tively for lensed quasars whose redshift of the lens galaxy is 
known and unknown) are consistent with the the expected color 
differences within the Icr error bar. This indicates that the ob- 
served chromaticity can be explained by differential dust extinc- 
tion. However, some data points cannot be reproduced by the 
currently applied extinction laws. For example, an upper right 
data point in panel (c) has error bars small enough to conclude 
that it deviates significantly from the theoretical lines. We will 
further discuss such data points in section 6. 

Here, we have applied only well-known empirical formulae 
for the extinction laws, and have not assumed any special, eccen- 
tric, or other unrealistic extinction laws. If it is confirmed that the 
differential dust extinction is responsible for the chromaticity of 
the lensed quasars, and if our assumption that the shape of ex- 
tinction curv e is the same at different positions in the lens galaxy, 
is valid (see, 'McGoug h et al.L l2005l) . we can conclude that any 
special extinction laws are not required to reproduce dust extinc- 
tion in galaxies up to z ~ 1, and the normal extinction laws that 
we have used in this study, i.e., the extinction laws derived from 
the local Universe, can be applicable in such distant galaxies. 



where /Imax and /Imin represent the maximum and minimum 
wavelengths of the filter response, respectively, and mi„,^k is the 
source magnitude without dust extinction. Finally, the expected 
magnitude difference between multiple images at filter k, Am^, 
due to differential dust extinction is expressed as 

Anik = mt(A^H) - '^^(A^h)' (5) 

where A^h and A^^ are the gas column density on one image and 
that on the other image, and the difference of these two quanti- 
ties is AA^H- Consequently, the expected color difference between 
multiple images is estimated by subtracting the magnitude dif- 
ference of one filter from that of the other filter, i.e., Am^ - Anij. 
Although the spectral shape of the source is clearly involved in 
equation m dependence on the assumed spectral shape causes a 
difference of only a few percent Q. Thus, we assume a flat spec- 
trum, i.e., S{A) - const. Applying equation |4] to all the filters, 

' By assuming a power-law spectral shape of the source, 5 (A) oc A^", 
we have checked the dependence. The expected dust extinction changes 
only a few percent by changing the power index (a) from to -3 



5. Quasar Microlensing 

From macro lensing modeling, the surface mass density of the 
lens divided by t he crit ical surface mass density for lensing (e.g., 
ISchneider et al.L 1 19921) is of order unity for all multiple images 
of lensed quasars. This indicates that the probability of quasar 
microlensing is high enough to be observed if the mass is dom- 
inated by compact components suc h as stars, planets or other 
compa ct objects. As investigated bv IWambsganss & Paczviiskil 
(Il991h . quasar microlensing can change the observed color of 
quasars from their intrinsic one, and yu(f) in equation [T] is also 
a function of wavelength. This effect has been examined by 
using more reahstic models for quasar accretion disks (e.g., 
lYoneharaetaTl 119981) . Since the spatial distribution of lens ob- 
jects on different images is generally different, quasar microlens- 
ing occurs in different ways in different images and consequently 
the color differences may be produced by quasar microlensing. 

This value could be too much as a representative value for the dif- 
ferential dust extinction, but the value indicates the maximum capability 
of the differential dust extinction to produce color differences. 
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Possible microlensing signals in the observed quasar sample are 
summarized in appendix B. 

Furthermore, the timescale of quasar microlensing can be es- 
timated by a timescale on which the source crosses the Einstein 
ring radius (fg) or a timescale on which the caustic crosses the 
source (tc). The former timescale is evaluated by 



tE ^ — 



1/2 



1 MGMiens DuDol 
V, \ Dos 



14 



/ Vt r' / Miens 

llO^kms-i/ \1M0 



1/2 



DisDqi/Dq 
1 Gpc 



1/2 



(6) 
(7) 



where r^, y, and Miens represent the Einstein ring radius on 
the lens plane, the transverse velocity of the lens on the lens 
plane , and the mass of the lens, respectively (e.g., Ilrwin et aLl 
119891) . The latter timescale is evaluated by 



tc = 



\10%ms-i/ \ 



Mu 



(8) 



where r^, Vc and Mbh represent the Schwarzschild radius, the 
transverse velocity of the caustics on the source plane, and the 
mass of supermassive black hole in the center of the quasar, re- 
spectively Q Of course, there are several ambiguities in these 
estimates, but the expected timescale for quasar microlensing 
should be comparable to or longer than 1 year. This timescale 
is long enough for the color difference caused by microlensing 
to be observed as a static phenomenon, and the expected color 
difference due to quasar microlensing has to be a realistic candi- 
date to reproduce the observed chromaticity. 

To quantify the microlensing effect, we should treat the 
continuum source (quasar accretion disk) and the magnifica- 
tion properties of microlensing. For a model of the continuum 
source, we applied the so-ca lled "standard accretion disk model 
jShakura & Sunvaevl Il973h " as a central engine of quasars. 
Here, the inner and the outer radii of the accretion disk are set 
to be 3r.5 and lO^r, 0, respectively, and the accretion rate is 
set to be a critical value. Thus, the only parameter specifying 
the properties of an accretion disk is the mass of the central 
supermassive black hole (Mbh)- Some emission Unes can con- 
tribute to the optical flux, and such emission lines can be af- 
fected by micr olensing in s ome case s (e.g ., Abaias et al., 2002; 
iLewis & Ibatal .2004: Richards et all 12004; Sluse et al.. 2007) . 
However, it is hard to imagine that all the sample quasars with 



' The transverse velocity of the caustics is determined by a combi- 
nation of a bulk and a proper motion of the lens objects. However, a 
proper motion of the lens objects dramatically changes the caustics net- 
works themselves, and the expected transverse velocity is lar ger than 
the simple combination of these two components (e.g., Wvit he et al.L 
11999 ). Moreover, how the caustics networks changes also depends on 
surface density of the lens objects, that of smooth matter, and external 
shear. 

The effective temperature of the accretion disk at a radius of 10' r, 
is ~ 2600 K, and the peak wavelength of black body spectrum with this 
temperature is ~ 1.1 x 10"* A. Assuming the source redshift to be Zs = 
1 .0, the peak wavelength at observer is ~ l.lxlO''x(l+Zj) = 2.2x10'* A. 
Since this is longer than wavelength coverage of the reddest filter in our 
calculations (F160W of HST NICMOS), the outer radius larger than 
~ 10' Tj does not change the results. If the radius is much smaller than 
lO'r,, the resulting microlensing signal tends to be enhanced. 



various redshifts are significantly affected by line emissions in 
the observational bands. Here, we neglect such line emissions 
and consider only continuum emissions. For magnification prop- 
erties, we have considered the properties of the caustics, which 
coiTesponds to regions where the source is extremely magnified 
by microlensing. In this paper, we take a straight line approxima- 
tion for fold caustics and apply an approximated magnification 
in the vicinity of caustics w hich is only a fun ction of a distance 
from the caustics (e.g., Sch neider et all [T992I) . Including a con- 
stant magnification, (magnification due to all the lensing effects 
except quasar microlensing), the magnification applied in this 
study is expressed as 

/ X Y"^ 

fj.(x) = — + fia (inside , or positive side of, the causti(9) 
- fia (outside , or negative side of , the caustics), (10) 

where x, Xs, and fj.a represent the distance from the caustics, the 
scale length of the caustics and the constant magnification, re- 
spectively Since this formula is investigated from general 
(mathematical) properties of gravitational lensing, this approx- 
imation is applicable not only for caustics produced by single 
lens object, but also for more general case such as caustics of 
quasar microlensing, i.e., caustics produced by multiple lens ob- 
jects with external convergence and shear. Effects of the distri- 
bution of lens objects, external convergence and shear on magni- 
fication are hidden in the scale length, x^, in this approximation 
(see more details in Appendix C). 

Magnification due to macrolensing is different at different 
images, and /jg is also different at different images. Then, the 
expected magnitude at filter k at epoch f, ink, including quasar 
microlensing magnification is expressed as 



mk(t,fla) = 

mint,k - 2.5 log 



£r ■f''^'^^ Xi-ea 1^^^^ ~ Xc(Olmin)Q/(l+z,,)(x) dx 



where q(x) and Xc(0 represent emissivity distribution of the 
accretion disk at wavelength A and the 2-dimensional location 
of caustics on the sky at an observational epoch, respectively. 
Finally, the expected magnitude difference between multiple im- 
ages at filter k, Ami, due to quasar microlensing is expressed as 



Amk = mk(t,ida) - mk(t',ij'J, 



(12) 



where t and t' are an epoch on microlensing at one image and 
that at the other image, respectively, and they should be mea- 
sured at the rest-frame of lens galaxy. Again, the expected color 
difference between multiple images are estimated by subtracting 
the magnitude difference at one filter from the other filter, i.e., 
Anik - Amj. fia and yu', are a macrolensing magnification factor 
of one image and that of another image, respectively. The inner 
integral should be performed all over the source, and the term 
|x - Xc(r)lmin, minimum value for |x - Xc{t)\, corresponds to the 
distance to the caustics. In this estimate, like the estimation for 
differential dust extinction, we again take into account the filter 



" When the source crosses fold caustics and comes 'inside' the caus- 
tics, two images appear at corresponding critical curves. If the source 
moves to the other way around, i.e., the source goes 'outside' the caus- 
tics, the pair of images will disappear at the critical curves. Even a part 
of the caustics is approximated as a straight line, such property should 
remain unchanged and is involved in equation [TO] The first term in the 
right hand side of equation [9] represents magnification for the appear- 
ing/disappearing pair of images at the critical curves. 
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Fig. 6. In the lower panel, the expected quasar microlensing light 
curves at F555W, F8 14W, and F160W are presented by the solid, 
dotted, and dashed line, respectively. Abscissa is time steps in 
unit of r jvc, and ordinate is magnification factor in linear scale. 
At the beginning of these light curves, the source is completely 
outside the assumed caustics and the location of the source cen- 
ter respective to the caustics is -1.5 x lO^r,. At the end of these 
light curves, the source is completely inside the assumed caustics 
and the location of the source center respective to the caustics is 
+3.5 X IQ^Ts- In the upper panel, corresponding color changes (in 
mag) during the events are also presented by the solid (F555W 
and F160W) and dotted (F814W and F160W) Hne. z,, Mbh and 
Ha are set to be 2, IO^Mq, and 1, respectively. The Einstein ring 
radius for zi = 0.5 and for Miens - ^Mq at the source plane is 
applied for Xg. 



responses. An example of microlensing light curves for all the 
three filters are presented in Figure |6] Since there is no correla- 
tion among location of caustics on different images, Xc must be 
different at different images. 

Based on the above definitions and assumptions, we have 
calculated microlensing light curves in all the three filters for 
various parameters. The light curves show a clear waveband 
dependence in their shapes (see Figure |6l); at shorter wave- 
length, the expected magnification changes more dramatically 
and the maximum magnification becomes larger This is because 
the emission at a shorter wavelength comes only from a rela- 
tively inner compact region of the accretion disk compared to 
the em ission at a longer wavelength (e.g., Shakura & Sunvaev, 
[T973h . A smaller source is magnified more in a microlens- 
ing event because a large fraction of the source can locate in- 
side a strongly magnified region. Further, as shown in equa- 
tion [8] the timescale of microlensing event depends on the 
source size, and the ev ent observed at shorter wavelength is rapid 
(IWambsganss & PaczyhskiLII991l:lYonehara et aliri998l) . For all 
calculations of the light curves, the source center is located at 
-1.5 X 10^ from the caustics at the beginning and is located at 
+3.5 X lO^rj from the caustics at the end. Since we assume 10^ rs 
as the accretion disk size, the source is completely outside the 
caustics at the beginning, and is completely inside the caustics 
at the end. 



The range completely covers the most interesting epoch for 
quasar microlensing events associated with caustic crossing, and 
the coverage is enough for our current purpose. Subsequently, 
we pick up all possible combinations of two epochs in the light 
curve. This mimics an observation of two images in multiple 
quasars with different time delay. Subtracting the magnitudes of 
the brighter image at F160W from that of the fainter one in all fil- 
ters, and calculating the color difference between the images, we 
can estimate how much color difference between images is ex- 
pected from quasar microlensing. The expected color difference 
due to quasar microlensing in our calculations is determined by 
5 parameters; Mbh, Zs, Ha for two images to be compared (de- 
noted as jiA and hb), and Xj. Based on realistic estimations, Xs 
takes values around and we apply as a standard value of Xj 
(see Appendix C for more details about x,). Since is a func- 
tion of Miens, zi and Zj, physical value of Xs is determined by 
Miens, zi and Zs- The expected area of the color differences are 
shown in Figure |7] 

It is apparent from Figure [T] that the maximum color differ- 
ence depends strongly on Mbh, ip^A,I^B) and x,, and weakly on 
Zs- Within a reasonable range for Zs (e-g-, see FigurelTJ, the rest 
frame wavelength which corresponds to the observational wave- 
band varies by a factor of few. Thus, emissivity distributions for 
the observational wavebands at the different source redshift dif- 
fers slightly. Not only due to the smaller source size, but also 
due to the bluer spectrum of the accretion, the expected color 
change for smaller Mbh shows different features than that for 
larger Mbh, as shown by the dotted line in Figure |7] (a). For 
larger ^a, i-S-, larger yu^, larger hb or both in these calculations, 
the color difference is suppressed as shown by the dotted line 
in Figure |2](b). This can be understood from equations |9] [10] 
and[TT] Magnitude changes due to quasar microlensing are de- 
termined by a ratio of a magnification for inside the caustics 
(equation |9j to a magnification for outside the caustics (equa- 
tion [TOll rather than a magnification for inside the caustics alone 
(see equation fTTTi. Since the contribution of the term including 
X in equation |9] to the magnitude change is relatively smaller 
in the larger Ha case, the change due to microlensing becomes 
small, and consequently, the color change is smaller compared 
to the smaller Ha case. This means that the magnification due to 
macrolensing could be crucial to quantify how much color dif- 
ference would be expected from quasar microlensing. Of course, 
another major factor to estimate the color difference is a prop- 
erty of magnification pattern as shown with the dashed lines in 
Figure|2](b). If we apply smaller x^, only a part of the source is 
magnified by microlensing, and a smaller amount of color dif- 
ference is expected (see equation [TOli. 

Since the 'standard' parameter set in Figures [T] (a) and (b) 
can be a representative parameter set for actual quasar mi- 
crolensing events, the expected color is also plotted with obser- 
vational data in Figures |7] (c) and (d) for comparison. Quasar 
microlensing alone can also produce color differences up to 
~ 2 mag which is similar to the scatter of observational data. 
Further, the observed color difference is well reproduced by 
quasar microlensing within the error bars except some individual 
data points, for instance the upper right data point in Figure|2](c) 
and the lower left data point in in Figure |2](d). 

As is already noted in the case of quasar variabilities (sec- 
tion 3), the date of observatio ns is different in different wave- 
bands (e.g.. lLehar et al.L[2000l) . and we should take into account 
the interval among the multi-waveband observations for one ob- 
ject, Af. However, conversion of the timescale of microlens- 
ing used in numerical calculation, e.g., see equation [8] into an 
actual timescale involves some ambiguities. Here, we applied 



10 



Atsunori Yonehara et al.: Origin of Chromatic Features in Multiple Quasars 



|Af| < SOrJvc and |Af| < lOOrJvc as the time difference with 
respect to the date of observation at F555W and calculate the 
expected color difference between images. That is, denoting the 
observational epoch at F555W in the microlensing light curve 
for two images by Ia and tg, the observational epoch for two im- 
ages at F814W and F160W is + At and Ib + At, and + At' 
and tB + At', respectively, where At and At' are uncorrelated and 
all possible combinations for At and At' within the applied time 
difference are chosen for calculations. The maximum time dif- 
ference of |Af| < SQrJvc and |Af| < IQQrJvc corresponds to 
~ 1 yr and ~ 2 yr, respectively, in the case of - 10^ km s"' 
and Mbh = IO^Mq. The results are presented in Figure [8] As 
easily seen in Figure [8j the expected color difference with non- 
zero time difference among different wavebands (the dotted and 
the dashed lines in Figure [8^) is stretched toward an orthogonal 
direction to the major axis of the expected color difference cal- 
culated without time difference (the solid line in Figure[8^). An 
effect of time difference largely extends the theoretically allowed 
region of the expected color difference produced by quasar mi- 
crolensing. Consequently, most of the observed color difference 
of multiple quasars including objects whose observed color dif- 
ference cannot be reproduced by differential dust extinction are 
also nicely reproduced (Figure [51 (b))PI. 

Since we do not include any statistical properties of caus- 
tics such as distribution of the scale length of the caustics or 
clustering of the caustics the expected color differences can not 
reproduce distribution of the data points well Q On one hand, 
therefore, we cannot directly compare to the results for quasar 
microlensing to the distribution of observational data (e.g., data 
points shown in Figure [8] (b)). The expected color differences 
presented here (e.g., a closed curve shown in Figure [8] (b)) are 
the maximum color differences which can be produced within 
the framework of our calculations with a given Xs- 

On the other hand, the applied approximation for magnifi- 
cation due to quasar microlensing is surprisingly sufficient to 
reproduce the observed color differences. Our treatment, an ap- 
proximation for magnification based on general properties of 
gravitational lensing, is rather simple compared to more realis- 
tic methods such as ray-shooting, and does not include any fine 
structure of the caustics such as curvature of fold caustics, cusp 
caustics[3, or an overlap of caustics. We suspect that our simple 
treatment is enough to reproduce the observed color differences 
since the observed color differences are mainly produced by iso- 
lated caustics and/or caustics with larger size compared to fine 
structures of caustics. Even if caustics with a typical size smaller 
than the source size overlap with caustics with a larger size, the 
former (smaller) caustics add only a tiny fluctuation of magni- 
fication, and cannot produce significant amount of color differ- 
ences. Thus, we can interpret that our simple approach extract 



The expected color differences presented here are the maximum 
color differences which can be produced within the framework of our 
calculations. Since we do not include any statistical properties of caus- 
tics such as distribution of the scale length of the caustics or clustering 
of the caustics, the expected color differences can not reproduce distri- 
bution of the data points well. This could be the reason why the data 
points are clustering rather narrow region in Figure [8] (b) compared to 
the expected color differences. 

This could be the reason why the data points are clustering rather 
narrow region in Figure [8] (b) compared to the expected color differ- 
ences. 

Since cusp caustics require more special conditions than fold caus- 
tics (Schneider et al. 1992), fold caustics crossing occurs much more 
frequently compared to cusp caustics crossing. Thus, ignoring cusp 
caustics crossing only misses some rare events. 



only the caustics which contribute to the strongest magnifica- 
tion at a certain moment. This putative explanation should be 
checked by more sophisticated calculations in future. 

Even with quasar microlensing, it is clear that the observed 
color differences of some extreme objects which are located at a 
top-right corner and a bottom-left corner in Figure [8] (b) are not 
well reproduced. We will further discuss this issue in section 6. 

6. Discussion 

In this paper, we have examined the following three possibili- 
ties for explaining the observed chromaticity between multiple 
images of lensed quasars: (i) intrinsic variabilities of the lensed 
quasars, (ii) differential dust extinction in the lens galaxy, and 
(iii) quasar microlensing in the lens galaxy. 

6. 1 . Overall trends 

First, we have examined the intrinsic variabilities of the lensed 
quasars. Quasars generally have intrinsic variabilities, and there 
exists a time delay between multiple images of lensed quasars. 
Thus, the intrinsic variabilities of quasars could explain the ob- 
served color differences. However, we find that the expected 
value of the color difference is roughly one order of magnitude 
smaller than the observed color differences. Therefore, we can 
exclude the possibility that intrinsic variabilities are responsible 
for the observed color differences in the sample. 

As shown in Figures |5] and |2l both differential dust extinc- 
tion and quasar microlensing instead can reproduce the observed 
color differences. The color differences due to these two scenar- 
ios similarly reproduce the observational data points in Figures|5] 
and |2] However, once we take into account the time difference 
among the observations at different wavebands that usually hap- 
pens in an actual situation, the area covered by the expected 
color difference of quasar microlensing expands, and we can re- 
produce the observed color difference by quasar microlensing 
alone, as shown in Figure [8] These results are also summarized 
in Figure |9j objects whose observed color differences are not 
reproduced by a single scenario alone, either differential dust 
extinction or quasar microlensing, are plotted with the expected 
color difference by our estimations. Although there is an ambi- 
guity in converting a timescale used in our numerical calcula- 
tions into an actual timescale, quasar microlensing proves to be 
a better scenario to reproduce observed color differences com- 
pared to other scenarios, if only one scenario is responsible for 
the observed color differences. In an actual situation, it is more 
probable that both effects play a role at the same time in produc- 
ing the observed color differences. Combining these two scenar- 
ios (differential dust extinction and quasar microlensing), all the 
observed color differences presented in this paper seem to be 
nicely reproduced (see Figure |9]i; if we convolve the shaded re- 
gion in Figures |9] (a) and (b), we may be able to reproduce the 
observed color differences (the gray crosses) shown in Figure |9] 
which cannot be reproduced by either differential dust extinction 
or quasar microlensing. 

Unfortunately, it is not possible to discriminate which one 
is dominating the observed color differences with only one 
photometric observation as used in this paper. Since differen- 
tial extinction is a static phenomenon and quasar microlensing 
is a temporal phenomenon whose timescale is given by equa- 
tions [T] and/or |8] multiple or monitoring observations of lensed 
quasars at multiple wavebands are the only way to distinguish 
between differential dust extinction and quasar microlensing. 
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Multi-wavelength data taken at the same time or during negligi- 
bly short spans of observations compared to timescale of quasar 
microlensing are also required to suppress the color diversity 
caused by the time difference among the observations at different 
wavebands. 

Of special interest are the properties of objects whose ob- 
served color differences are not reproduced by either differen- 
tial dust extinction or quasar microlensing. Even if the observed 
color difference is within the expected color difference of our 
estimations, we have to care about quasar microlensing and/or 
intrinsic variabilities of quasars. Thus one should trace and 
separate effects of quasar microlensing and/or intrinsic quasar 
variabilities properly to derive the extinction curve of the lens 
galaxy: otherwise the extinction curves derived for lens galax- 
ies may be contaminated by the microlensing and/or the intrin- 
sic quasar variability effects, and it will be an unrealistic at the 
worst. 

6.2. Peculiar features 

Next, we focus on objects whose observed color differ- 
ences are difficult to be reproduced by either differen- 
tial dust extinction or quasar microlensing. In Figures |9] 
(a) and (b), we choose 8 objects (10 image pairs) whose 
color differences are not explained by differential dust ex- 
tinction alone (B0218H-357, SBS0909-H523, LBQS 1009-0252, 
PG1115H-080, SBS1520-H530, B2045-H265, Q2237-H0305, and 
APM08279H-5255), and 3 objects (3 image pairs) whose color 
differences are not explained by quasar microlensing alone 
(B0218+357, MG2016+112, and APM08279-H5255). 

A possible explanation for eight objects deviating from the 
expected color difference in the differential dust extinction sce- 
nario is that the lens galaxies have a peculiar dust extinction 
which cannot be parameterized by the function that we have 
used in this paper Since galaxies at different evolutionary stages 
are e xpected to have dust with different extinc tion properties 
fe.g.. lMaiolino et al.Ll2004HHirashita et al.Ll2005h. this wo uld be 
one plausible explanation. Furthermore, Inoue et aTl (12006 ) have 
shown that the wavelength dependence of the dust attenuation 
is modified by a effect of light scattering. An anomalous dust 
extinction (or attenuation) might be seen in a systematic differ- 
ence in the colors of lens galaxies themselves. Thus, in order to 
check if the anomalous dust extinction is really responsible for 
the deviation from the model predictions, we examine the colors 
of the lens galaxies in Figure [TO](see also Table|2]i. Photometric 
data for the lens galaxies of B1600+434 and Q1208-H101 are 
available not for all the bands, and the data are not plotted 
in Figure[TO] B1030-H071, MG2016-I-112, FBQ095 1-1-2635, and 
Q0957 -1-561 have another lens candidate, and the data are also 
plotted in Figure [TOl 

Unfortunately, as we can see in Figure [TOf a). this explana- 
tion does not seem to be confirmed. Colors of the lens galaxy of 
objects with the unexpected color differences and that with the 
expected color differences are compared in Figure fTOl a). The av- 
erage color for the former objects are {mfs^sw-'nFmw, mpsuw- 
mpieow) - (3.93 + 1. 12, 2. ll+0.74),and that for the latter objects 
are (otf555w - '"Fieow, '«f814w - '«fi60w) = (4.35 + 1.40,2.33 + 
0.53). The average color of the lens galaxy of objects with the 
unexpected color differences is slightly bluer than the others, 
but these values are overlapping each other within the disper- 
sion. Thus, we can say that these two groups are not separated 
clearly as far as the galaxy color is concerned. Further, there is 
no clear difference between the redshifts of these two groups. 
The average redshifts are (zi, Zs) = (0.62 + 0.30, 1 .94 + 0.88) and 



(zi,z,) = (0.61 + 0.23,2.36 + 0.97) for the objects with the un- 
expected color differences and the others, respectively (see also 
Figure[TO]b). 

In principle, dust extinction should be larger at shorter wave- 
length. This limit corresponds to the dashed lines in Figure |9l 
the horizontal line indicates that absorption at F814W and that 
at F160W is equal, and the diagonal line shows that absorption 
at F555 W and that at F8 14W are equal. Therefore, the color dif- 
ferences by dust extinction should always be between the dashed 
lines in Figure |9j a). In other words, the color difference of ob- 
jects that are not located in this region is not reproduced by dust 
extinction alone. As shown in Table |2] some of the lens galaxies 
are known to have observational properties of late-type galax- 
ies; B0218+357, RXJ0911-H0551, B1030+071, SBS1520+530, 
B1600H-434, B2045-H265, and MG2016-I-1 12, and dust proper- 
ties in such systems may be somehow different from that in the 
MW or the SMC. It is still possible to explain the color differ- 
ences in B2045 -1-265 by differential extinction of dust with spe- 
cial properties because the color diflFerences are located at "al- 
lowed" area for differential dust extinction. However, color dif- 
ferences in B0218+357, SBS 1520-1-530, and Q2237+0305 are 
located at "prohibited" area for differential dust extinction, and 
their anomalous color differences are hard to be explained by 
differential dust extinction alone, if dust extinction is assumed 
to cause reddening in any wavelength. However, as can be seen 
from Figure [TO] there is no clear evidence for peculiarity in the 
colors and redshifts of the objects whose deviation from the 
models is large with respect to other objects. 

As for the 3 objects deviated from the expected color dif- 
ference in the quasar microlensing scenario, it is possible that 
the lensed quasar have a peculiar physical properties of an ac- 
cretion disk which cannot be well represented by the standard 
accretion disk model that we have used in this paper Since the 
expected color differences by quasar microlensing may depend 
on the radiative properties of the central engine of quasars, a sys- 
tematic difference in the source color between these two groups 
would be expected. In Figure [TOf c). we present the color-color 
diagram of all the source quasar images. Again, however, no 
clear difference between these two groups is found. The aver- 
age colors of the source quasars of the systems with unexpected 
and expected color differences from our microlensing models are 
(mF555w -mFi60w,mFm4w -niFieow) - (3.76 + 1.08,2.37 + 1.00) 
and {mF555w -mFmw,mFm4w -mFi60w) - (3.03+2.31,1.75 + 
1.16), respectively. Thus, as long as we focus on the colors of 
the source quasars, we cannot find any significant difference 
between these two groups. Further, there is no clear difference 
between the redshifts of these two groups. The average red- 
shifts are (z,,z,) = (0.85 ± 0.17,2.70 + 1.25) and (zi,Zs) = 
(0.58 ± 0.24,2.16 + 0.90) for the objects with the unexpected 
color differences and the others, respectively (see also Figure [TOl 
d). 

A part of objects whose observed color differences are dif- 
ficult to reproduce has been observed with more wavebands 
than others. For instance, B02 18+357 and APM08279+5255 
have been observed with 6 and 5 different wavebands, respec- 
tively. By using those data, we can examine the color dif- 
ferences between images more intensively. In Figure [TT] we 
present the magnitude differences between images. In the case 
of APM08279+5255, the trend is rather simple. The magnitude 
difference monotonically decreases with increasing wavelength, 
and bluer photons are absorbed more than redder photons at im- 
age B and/or are magnified more than redder photons at image 
A. As we can also see from Figure |9] the color difference can 
be reproduced by a flatter extinction law and/or a larger scale 
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length of caustics than that we used in this paper. The extinction 
law becomes flatter if the grain size distribution is more biased 
to a larger size. Furthermore, a larger scale length of caustics 
is also realized by taking into account large variety of caustics 
(e.g., Appendix C), and thus, it is possible to reproduce observed 
color differences of APM08279+5255 either by dust extinction 
or by quasar microlensing. However, the trend is complicated in 
the case of B0218+357, because the observed color difference 
seems to have a peak around 7000A. If we can treat two data sets 
presented in Figure [TT]equivalently, the magnitude difference is 
0.5-1.0 mag at ~ 3000A and at ~ 16000A, but is ~ 2.5mag 
at ~ 7000A. It is clear that the magnitude difference rapidly in- 
creases up to a wavelength of ~ 5000A and becomes almost con- 
stant at wavelengths of ~ 5000 8000A. Such a wavelength 

dependence cannot be reproduced by dust grains. In contrast, 
quasar microlensing can magnify a part of quasar accretion disk 
selectively, and an emitting region of photons with a wavelength 
of ~ 7000A can be selectively magnified with a certain con- 
figuration of the accretion disk and caustics. Further, if photons 
with wavelength above ~ 7000A mainly originate from more 
extended region than the accretion disk, such as a lobe of jets, 
a torus or a starburst region, these photons cannot strongly be 
magnified by microlensing. Considering this possibility of mi- 
crolensing, B0218-t-357 is more likely to be magnified by quasar 
microlensing than dust extinction. 

As for the observational properties of objects whose ob- 
served color differences are not reproduced by either differential 
dust extinction or quasar microlensing, we cannot find any clear 
difference between these objects and the objects whose observed 
color differences are well reproduced by either differential dust 
extinction or quasar microlensing. Therefore, the peculiarity of 
lens galaxies or lensed quasars is unlikely to be a probable reason 
for the extreme color differences. It may be more natural to con- 
sider that all the observed color differences can be reproduced 
by a combination of all three possibilities presented in this paper 
without requiring any special properties of dust and/or quasar 

6.3. Concluding remarks 

It is worth mentioning that there are some limitations in our 
current treatments of differential dust extinction and quasar mi- 
crolensing. 

For differential dust extinction, we have applied empirical 
extinction laws derived for the local universe such as the MW 
and the SMC. Although the observational color differences are 
roughly consistent with the local extinction curves, the extinc- 
tion laws in distant galaxies such as the lens galaxies may be sig- 
nificantly different from the local ones, especially in lens galax- 
ies whose data points deviate significantly from our calculations 
in Figure |5] In contrast, we have used gas distri bution obtained 
by a numerical simulation for a disk galaxy (see lHirashita et all 
12003 ) to investigate the expected differential dust extinction. The 
current sample of lens galaxies involves both early- and late-type 
galaxies (see Table |2]l, which may have very different gas and 
dust properties. To take into account such variety of lens galax- 
ies, gas and dust distribution for various types of galaxies, espe- 
cially for early-type galaxies, and further exploration by using 
the distribution is required. 

For quasar microlensing, we applied a straight line caus- 
tic approximation and did not take into account com- 
plex caustic networks ( e.g., Wambsganss & Paczvriski, 119911 : 
iGoicoechea. et al.ll2005h that would be expected in most of sys- 
tems. On the one hand, multiple caustics will produce more dra- 
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Fig. 11. Magnitude difference at several wavebands between im- 
ages of B0218+357 (upper panel) and APM08279-t-5255 (lower 
panel). These data are also provided by CASTLES web page. 
For B0218-I-357, magnitude differences (magnitude of image 
A minus that of image B, the brightest image at F160W of 
HST NICMOS) at F300W, F380W, F439W, F555W, F675W, 
and F814W of HST WFPC2 are presented by the solid eiTor 
bars. Data used previously in this paper are different sequence 
from these photometric data at 6 wavebands, and the data used 
previously are also presented by the dashed error bars as refer- 
ences. For APM08279+5255, magnitude differences at F555W 
and F814W of //5rWFPC2, and at Fl lOW, F160W, and F205W 
of //5r NICMOS are presented. The solid error bars indicate that 
magnitude of image B minus that of image A, the brightest im- 
age at F160W. 



matic and/or complex structures in the expected light curves 
of quasar microlensing, and even larger or more strange fea- 
tures will be expected in the observed color difference. On the 
other hand, caustics with short scale length, i.e., small jCj, or 
small caustics can just slightly change the obser ved color during 
microlensing events. As already pointed out bv lYonehara et al.l 
(Il998h . the expected color difference also depends on accretion 
disk models for quasars, because the spatial emissivity distribu- 
tion and the wavelength dependence is diverse in different ac- 
cretion disk models such as a r adiatively inefficient a ccretion 
flow. For instance, as presented in lYonehara et al.l ( Il998h . the ex- 
pected light curves of microlensing of an advection-dominated 
accretion flow (ADAF) is almost achromatic, at least in optical 
wavebands, since the spatial emissivity distribution of ADAF 
is almost independent of the wavelength. Thus, the ADAF 
model has difficulty in explaining the achromaticity of lensed 
quasars. Though the dependence on the source emissivity pro- 
file is less prominent in statistical properties of the light curves 
(iMortonson et al.Ll2005b . there exists clear dependence not only 
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on the source size but also on the source emissivity profile in in- 
dividual light curves. Quantitative estimates of various complex- 
ity caused by a variety of extinction curve, of caustic networks, 
and of accretion disks are addressed in future works. 
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Fig. 5. The upper two panels show the expected color differences between multiple images for various extinction laws at the same 
lens redshift (panel a), and for the same extinction law at various lens redshifts (panel b). All of these lines are drawn up to the 
expected color differences for 95% of the cumulative differential column density of the gas in the lens galaxy, AA^h ~ 7.1x10^' cm"^. 
In the upper left panel (panel a), the expected color differences for the MW extinction curve with Rv - 2.6, with Rv = 3.1, and 
with Rv = 5.6, and for the SMC extinction curve are presented by the solid, dotted, dashed, and long-dashed lines, respectively. 
The applied redshift for the lens in these estimations is 0.5. In the upper right panel (panel b), the expected color difference in the 
case of the lens galaxy at zi - 0.0, 0.5, 1.0, and 1.5 are presented by the solid, dotted, dashed, and long-dashed lines, respectively. 
The applied extinction law in this estimation is the MW one with Ry = 3.1, or a standard extinction law in our galaxy. The lower 
two panels show the expected color differences for the various extinction laws at the lens redshift of zi - 0.5. The extinction law for 
the MW with Ry - 2.6, the MW with Ry - 5.6, and the SMC are presented by the solid, dotted, and dashed lines respectively. The 
observational data which were shown in Figure|2](a) and (b) is overlaid as the gray crosses in panels (c) and (d), respectively. Note 
that the scale of the abscissa and ordinate in the upper panels is different from that in the lower panels. 
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Fig. 7. The expected color differences between multiple images due to quasar microlensing. The maximum color differences by 
quasar microlensing are presented by the lines. The upper two panels show the expected color difference for various parameter 
sets; Mbh, Zs, I^a, fJ^B, and x,. The solid line in all panel represents the expected color difference for Mbh = IO^Mq, Zs = 2.0, 
- (3.0,2.0), and Xs = r^, and is referred to as 'standard' in all the panels. Miens - ^Mq and zi - 0.5 are applied for 
calculating in all the cases, and in this estimation is not on the lens plane, but on the source plane. In the upper left panel, 
panel (a), the expected color differences for Mbh = lO^M© and for Zs = 1-0 are presented by the dotted (red) and dashed (blue) 
lines, respectively. In the upper right panel, panel (b), the expected color difference for {ij.a,I-Ib) = (10,2) and for x, = O.lr^ are 
presented by the dotted (red) and the dashed (blue) line, respectively. The remaining parameters are the same as that for parameters 
of the 'standard' case. Observational data shown in Figure|2](a) and (b) is plotted by gray crosses in panels (c) and (d), respectively, 
and the expected color differences for Mbh = 10'^ Mq presented in panel (a) are overlaid on these panels by the solid line. 
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Fig. 8. The expected color differences between multiple images due to quasar microlensing. Here, possible ranges of observations are 
included. 'Standard' set of parameters in Figure |7] are also applied for these calculations, and the result obtained without including 
the time difference, i.e., |Af| = 0, is presented as a reference by the solid line in the left panel, panel (a). In panel (a), the result for 
|Af| < SQrJvc and |Ar| < lOOrj/vc is also presented by the dotted and the dashed line, respectively. Observational data shown in 
Figure|2](a) and (b) are plotted by gray crosses in panels (b). The expected color difference in the case of |Af| < IQQrJvc is also 
overlaid on this panel by the sohd line. 
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Fig. 9. Objects whose observed color difference between multiple images is not well reproduced by our estimations are picked out 
and are presented with their name, and with the expected color difference by our estimations. In the left panel (panel a), the expected 
color difference for the MW extinction curve with Ry - 2.6 and that for MW with Ry - 5.6 are presented by the solid and dotted 
lines, respectively, where zi = 0.5 is assumed. The expected color difference locates between the solid and the dotted line, i.e., 
the shaded region. Objects whose observed color difference cannot be reproduced by differential dust extinction are also presented 
by the gray crosses. In the right panel (panel b), the expected color difference for 'standard' case of quasar microlensing and that 
including difference in observational epochs at different wavebands are presented by the solid and dotted lines, respectively, where 
Zs = 2.0 is assumed. The expected color difference for 'standard' case with the non-zero spans locates inside the solid line, i.e., the 
shaded region. Again, objects whose observed color difference cannot be reproduced by quasar microlensing are presented by the 
gray crosses. The line for Aiufkhw - Awf leow = and the line for AmF555w - AniFieow - Atofsmw - AniFieow are also presented 
by the dashed lines as references. 
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Fig. 10. Relations between the quantities concerning the lens galaxies and the source quasars. The black filled squares in the upper 
and lower panels represent the objects that have deviation from the expected color difference by differential dust extinction and by 
quasar microlensing, respectively. Properties for the remaining objects are also plotted by the gray filled triangles. The color-color 
diagram for the lens galaxies and the source quasars are presented in panels (a) and (c), respectively. In panels (b) and (d), the lens 
redshift and the source redshift are presented in the same way as panels (a) and (c). The objects with unknown lens redshifts are 
plotted on the diagonal dotted line. 
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Appendix A: Extinction Law in the SIVIC 

iGordon et alj (l2003l) have investigated extinction curves of the Small Magellanic Cloud (SMC), the Large Magellanic Cloud, and the 
Milky Way by using photometric data in the wavelength range from ultraviolet to near-infrared. Since the coefficient Ci in equation 
(5) of^ordon et al. ( 2003) is negative (Table 3 in Gordon et al., 2003) for the SMC, the extinction at some long wavelength becomes 
negative value, i.e., A{x ^ 0) < for positive Ay- This may not be ideal in some situation, and we ma de a fitting formulae for it by 
ourselves by using the same data for the SMC (denoted "SMC Bar" in Table 4 of lGordon et al.Ll2003l) . 

We have adopted the same functional form of the extinction cur ve as the one provide d by equation 1 in lCardelli et alj ( Il989h in 
our fitting procedures. We apply the same function form as a(x) in ICardelli et al.l (Il989h . and the fit is satisfactory as shown later 
The number of parameters to be determined is 20 in total; 1 coefficient and 1 power index at /I ' < 1.1 /im 7 coefficients at 
1.1 /im"' < /I ' < 3.3 /im"', 7 coefficients at 3.3 /im ' < /I ' < 8.0 ;um~', and 4 coeffi cients at 8.0 /zm ' < /I We have determined 
those coefficients by performing minimization with a downhill simplex method (P ress et al.Lll992l) . 

The resultant formulae for the extinction as a function of the inverse of wavelength (x = in yum"^) is as follows; 



A(x) 
Av 



= 0.31884x"'^'" (x < 1.1) (A.l) 

= 1 + 0.92733/9 - 0.31777/ - 0.83093/?^ + 1.58876/ 

+0.60622/- 1.91599/ + 0.71769/ (1.1 < x < 3.3) (A.2) 

0.02248 



-0.77359 + 0.8 1739X + 



(jc- 4.54152)2 +0.13635 



+0.23260/ - 0.06653/ (3.3 < x < 8.0) (A.3) 
= 6.17700 + 2.39744r- 8.90103r2 + 11. 44564r^ (8.0 < x), (A.4) 

where p - x - 1.82, q - x - 5.9, and r - x - 8.0. The;i'2 value for the best fit parameters is 36.28. Since the number of data points 
that we have used in our fitting is 30 and the degree of freedom becomes 30 - 20 = 10, the reduced ;t'2 value clearly exceeds 1 .0 and 
thus is not particularly good. This mainly originates in a few data points with small error bars at longer wavelength, and we can see 
in Figure |AT| that our fitting result nicely traces most of the data points, especially in the wavelength range relevant in this paper. 



Appendix B: Supplementary remarks for individual objects 

B. 1 . Possible evidence for quasar microlensing 

From previous measurements, the existence of a quasar microlensing signal is quite obvious from their observed light curves. For 
instance, Q2237+0305 has frequentl y shown cle ar microlensing signals (e.g.. IWozniak et all 12000) in the photometric data since 
the first detection of microlensing by llrwin et al.l (11989). Additionally, due to the large leverage arm of the nearby lens galaxy, the 
transverse velocity of the caustics on the source plane is enhanced by a factor of 10 from that of other lensed quasars, and the typical 
timescale of the microlensing events is reduced by the same factor. 

Po ssible microlensing signals in optical wavebands are sugg ested in Q0957+561 b y'Collev & Schild' ('2000'); Goicoechea. et al.l 
(120051) . and ones in radio wavebands ai'e suggested in B1600+434 by lKoopmans et al . (2001). Burud et al. (2000) have also shown 
a possibility of quasar microlensing signal in the /-band monitoring data of B 1600+434. 

For the measurement of time delay between images, contributions from quasar microlensing are included in the analysis, and the 
goodness of fit is improved substantially by including effects of quasar microlensing. F rom this point of view, the existe nce of quasar 
micro lensing is indicated indirectly in a subset of the cuiTent sample, PGl 1 15+080 (Barkanal[l997|), RXJ091 1+0551 (Hjorth et all 
120021) . HEl 104-1805 (Ofek & Maoz, 2003), SBS1520+530 (GaynulHna. et al., 2005), and FBQ095 1 +2635 dlakobsson et al., 2005F 
Observations of all these objects are performed at optical wavebands. 

For LBQS 1009-02 5 2, sparse multi-color m onitoring observations have been done and some spectra have been obtained by 
dClaeskens et al.Ll2001h . lClaeskens et al.l ( l200i1) have focused on the changes of colors and spectra during the observations and have 
mentioned a possible existence of dust extinction and/or microlensing in this system. However, it may be possible to explain the 
changes by intrinsic variabilities. 

In the case of QO 14 2- 100, evidence for quasar microlensing is considered to be obtained from the observed color changes in 
optical monitoring data dNakos et al.Ll2005h . 

Also in the case of HE2149-2745, quasar microlensing is suggested to explain unexpected spectral differences between images 
dBunid et aUl2002al) . 



B.2. Otiier remarks 

One of the most famous lensed quasars, Q0957+561, is in a somewhat special environmen t compared to other objects, because the 
lens galaxy is obviously a member of a galaxy cluster (e.g.. lAngonin-Willaime et an,ll994l) . 

MG2016+112 is supposed to have two lens galaxies. Although the lens galaxies show the brightness profile of early type 
galaxies, one or both of the lens galaxies show an [O ii] emission line (Koopmans & Treu, 200^ . 

Radio monitoring observations of MG0414+0534 have been performed bv lMoore & Hewitt! dl997l) . Unfortunately, the intrinsic 
flux variation of this quasar seems to be small, and a time delay measurement has not been successful. 
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Fig.A.l. The SMC extinction curve obtained by our fitting procedure is shown by the gray solid line. Abscissa is inverse of 
wavelength (/I '), and ordinate is extinct ion value divided by extinction at V band (A^/Ay). Data points with the error used in our 
fitting procedure (originally presented bv lGordon et al.Ll2003l) are overlaid as the black bars on the extinction curve. 



B 1030+074 has been detected bvl Xanthopoulos et al.l ( ll998 h. and at the same time, they claimed that an image of the lens galaxy 
which is taken by HST shows substructure which could be a spiral arm or an interacting galaxy. 

For B2045+265, all the three images show almost the same color difference from, or redder than, the brightest image in F160W 
filter (see Figure |9]l. However, it is not probable that three of four images are occasionally affected by dust extinction and/or 
quasar microlensing in the same way. Rather, it would be reasonable to consider that only the brightest image in F160W filter is 
predominantly affected by dust extinction and/or quasar microlensing, which cause a common color difference for the remaining 
three images. 



Appendix C: Scale Length of Caustics 

To estimate color changes produced by quasar microlensing, we have applied a simple approximation in the vicinity of fold caustics 
shown in equation|9l This approximation formula includes a scale length which is denoted by x,, and this scale length is not arbitral 
but determined by a property of the caustics (e.g., Schneider et al., 1992), i.e., external convergence and shear, and distribution of 
lens objects. Here, we briefly investigate how to estimate the length, and we present the expected value. 



C. 1 . Definition of tine scaie iengtin 

By using a lens potential, the lens equation for quasar microlensing with point mass lenses is expressed as 



= ^ = -A. + (1 - Kc)e, + r(cos 24,0, + sin 2(p0y) - 



N I r/] \2 



'del 



N 



a'' 



= -Py + (1 - Kc)0y + y(sin 200, - cos 200,,) - ^ 



|6/-6/['l| 



(0. - e?) 
(0, - Of), 



(C.l) 
(C.2) 
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where Px, Pv, Gx, Sy, ^i'' , ^l' ' , a^'^ , Key, and ^ represent x- and y- coordinates of source, x— and y— coordinates of image (0),x- and 
y- coordinates of i-th lens (0'''), Einstein ring radius of i-th lens, convergence due to smooth matter, external shear and direction of 
the shear, respectively. For a practical convenience in deriving the approximation formula, the coordinate systems are chosen such 
that the origin of the lens plane and the source plane is located on critical curves and caustics, respectively, and y-axis of the source 
plane is orthogonal to the fold caustics. Performing the Taylor expansion and applying conditions in the vicinity of the fold caustics, 
we can obtain the following expression for magnification (fj) after some algebra. 



1 


'5V 




2 „o 


A' 

(0,0) 


2 


[del 


(0,0). 


56f 



-1/2 



(C.3) 



Here, the caustics are assumed to be straight lines, and the dependence on yS, disappears from the magnification. This is the same 
form as the first term in equation |9] for inside the caustics, and now it is clear that the scale length in equation |9] should be 



X, - 2 



( 5V 








[ dOl 


(0,0). 


del. 


(0,0) _ 



(C.4) 



See lSchneider et al.1 (Il992h for more details. 



C.2. Efficient method to estimate fhe scale lengtli 

Altho ugh the magnification probability of quasar microlensing depends on the mass spectrum of the lens objects (ISchechter et all 
12004) . we assumed that all the lenses have the same mass. Further, all the length scales (jc,) are normalized to the Einstein ring 
radius for an assigned lens mass, and a''' is set to be unity here. The following estimation method can be applicable also if the mass 
spectrum is taken into account. 

In addition to equations lC.ll and lC.2l the following conditions should also be satisfied on fold caustics, or corresponding critical 
curves, i.e., at (0,, 0,,) - (0, 0); 

= ^-7- =rsm20+ > ■ (C.5) 

dxdy le-e"! 
= -^^(i-^,)-rcos20 + y 1^ — '-r — — — (C.6) 

^ -^ = (l-^,) + rcos20+ V . -' (C.7) 

dx^ jrt 

^ ^ = y 2(0,, - Of)- ^ —. (C.8) 

As we can recognize in the above equations, equation lC.7l is simply reduced to 2(1 - Kc) by using equation lC.6l Therefore, the only 
term we should evaluate to estimate Xs is the third derivative expressed in equation lC.81 A condition for equation IC. 8 1 corresponds 
to a condition to avoid cusp caustics. 

Since the above equations should be satisfied on fold caustics or corresponding critical curves, i.e., (6*1,6*,,) = (0,0), equa- 
tions |CT1|C21|C3] and lC.6l can be rewritten as the following forms. 







0[N-l]p 


|0[N]|2 


qIN-I] 


0^ 



/■=! 



N-2 9fl[/]fl[i] jdN-^AN-l] r.AN]alN] 

ZZUy Uy, ZU y Uy AU y Uy 

+ ^ + — ^ (c.ii) 

,=1 jeWf \0y^^\ 
= (1 -/fc)-rcos20+ ) — - — + — ■ + — . (C.12) 

Distributing the lens objects with the index from / = 1 to / = - 2 randomly, equations IC.9I IC. 1^ IC. 1 1 1 and IC.12l are equivalent 
with 4 equations for 4 unknown quantities (0^^"'^, 6''/'', and 6'[.^') for given y and (p. By solving those 4 equations, we can 

obtain locations of two lens objects to create fold caustics on the origin of the applied coordinate. Substituting locations of all the 
lens objects including 6'[^"'1 and 6''^' into equation |C.8l we can obtain for the fold caustics in this realization of distribution of the 
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Table C.l. The scale lengths of caustics, Xj, at certain cumulative probabilities for distributions presented in Figure lCTI are listed. 
The column for 50% indicates the median value of the distributions. 15.85% and 84.15% indicate the +lcr confidence levels from 
the median value, and 0.135% and 99.865% indicate the +3cr confidence levels from the median value. 
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lens objects. If the number of the lens objects which are randomly distributed on the lens plane (A^ - 2 in this case) is large enough 
compared to the number of the lens objects which are added to create the fold caustics at a certain place on the source plane (2 in 
this case), the resultant distribution of the lens objects can be also treated as a random distribution. 

An advantage of this procedure is that calculations of caustics networks are not required, and that computations for one realiza- 
tion are very fast. Although the procedure can only provide x., of caustics rather than the whole magnification pattern, it is efficient 
to investigate statistical properties of x,. The procedure must be useful to evaluate a proper value of Xj for a given environment, i.e., 
total convergence (k), Kc, and y. 

We have estimated x^ with 10"* lens objects (i.e., A^ = 10"*) for several values of k, k^, and y. All the lens objects except the last 
two, (A^- l)-th and A^-th lens objects, are randomly distributed within a circle. Radius of the circle is determined so that the surface 
density of the lens objects is equal to a- - Kc, i.e., [N/(k - /Cc)]'^^- Direction of external shear, (p, is also chosen randomly at every 
realization. We have performed 10^ realizations to obtain the expected distributions of x,, and results are presented in Figure ICTl 

First, the distribution depends on the fraction of smooth matter with respect to the total convergence, KcIk. The scale length x, 
becomes smaller as KcIk becomes smaller Since KcIk cannot be smaller than 0, the dashed line in Figure ICTT a) can be treated as a 
minimum of x, distribution for given k and y. Although it is practically difficult to investigate KcIk, we obtain the minimum of x^ for 
K - 0.5 and y - 0.5 as ~ 10-^ from Figure |C.ir a)PI. Second, as is shown in Figure ICTt b) and (c), dependence on k and y is weak. 
The minimum value of x., increases slightly with decreasing k, and the maximum value of x, increases slightly with decreasing y. 
Especially, the dependence of the Xj distribution on y is almost negligible. Although the distribution of x, shows a slight dependence 
on KcIk, k, and y, we can conclude that typical value of x^ is within a range of 0.1 - 1 as shown in Figure lCTI Several scale lengths 
at certain cumulative probabilities are also summarized in Table lC.ll for convenience. 



'^ The results may show a different aspect of what ISchechter et alj ( l2004h found. However, the distribution obtained here cannot be simply 
transformed into magnification probabilities, because we did not take into account neither spatial distribution of caustics nor size of caustics. 
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Fig. C.l. Cumulative probability distributions of the scale length of fold caustics, Xj. Abscissa is Xs in units of the Einstein ring 
radius (/"£■). In the top panel (panel a), total convergence (k) and external shear (y) are fixed to 0.5. Distributions for a fraction of 

smooth matter (kc/k) equals 0.0, 0.5, and 0.9 are presented by the dashed, solid, and dotted lines, respectively. In the middle panel 
(nanel h). k^I k and y are fixed as 0.5. Distributions for k = 0.8. 0.5. and 0.2 are presented bv the dashed, solid, and dotted line. 
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Table 1. Current sample of lensed quasars. The redshifts of the lens galaxy (z/) and of the source quasar (zs), the colors of the 
lensed quasar images and the errors for F555W - F160W (wfsssw - /WFieow). and for F814W - F160W (niFsuw - mprnw) are 
presented. These colors should be identical for each image of the same lensed quasar unless intrinsic variability, dust extinction, or 
microlensing play a significant role. The magnitudes of the images at F160W (m^ igoiv) are also presented in the last column. The 
redshifts of the lens galaxies are not known for the 6 objects listed at the bottom. All the data are taken from the CASTLES web 
page, and the colors are calculated from these data. 



Object Name 


Zi 


Zs 


r J J J vv 


" /71 F1 AOW 
r I uv) vv 












Q0142-104 A 


0.49 


2.72 


1.61 ± 


0.092 


1 


.34 


+ 


0, 


.181 


15 


.28 


+ 


0, 


.02 


B 






1.56 ± 


0.050 


1 


.14 


± 





.153 


17, 


.57 




0, 


.03 


B0218+357 A 


0.68 


0.96 


5.76 + 


0.036 


4 


.31 


+ 





.073 


17, 


.52 


+ 


0, 


.02 


B 






4.17 + 


0.022 


2 


.45 


+ 





.028 


16 


.94 


+ 


0, 


.02 


MG0414+0534 Al 


0.96 


2.64 


9.99 + 


1.270 


4 


.89 


+ 





.063 


15 


.54 


+ 


0, 


.02 


A2 






10.82 i 


: 0.272 


5 


.49 


± 


0, 


.067 


15, 


.87 




0, 


.03 


B 






9.52 ± 


2.080 


4 


.68 


± 





.076 


16 


.56 




0, 


.03 


c 






9.27 + 


0.141 


4 


.69 


+ 





.073 


17 


.41 


+ 


0, 


.02 


B07 12+472 A 


0.41 


1.34 


3.78 ± 


0.277 


2 


.64 


+ 





.434 


20, 


.46 


+ 


0, 


.21 


B 






4.58 + 


0.231 


3 


.14 


+ 





.422 


21 


.08 


+ 


0, 


.22 


c 






5.25 ± 


0.272 


3 


.32 


± 


0, 


.467 


21, 


.05 




0, 


.08 


D 






4.77 ± 


0.940 


2, 


.97 


± 





.433 


21 


.90 




0. 


.14 


RXJ091 1+0551 A 


0.77 


2.80 


1.24 + 


0.045 





.79 


+ 





.045 


17, 


.59 


+ 


0, 


.02 


B 






1.53 + 


0.054 





.99 


+ 





.036 


17 


.65 


+ 


0, 


.02 


c 






1.55 ± 


0.067 




.02 


+ 





.050 


18 


.34 


+ 


0, 


.03 


D 






1.62 ± 


0.102 


1 


.01 


± 





.036 


18 


.65 




0, 


.02 


SBS0909+523 A 


0.83 


1.38 


2.53 ± 


0.082 


1 


.95 


+ 





.028 


14 


.60 


+ 


0, 


.02 


B 






3.46 + 


0.173 


2 


.20 


+ 





.222 


14, 


.73 


+ 


0, 


.03 


BR10952-0115 A 


0.63 


4.50 


3.59 ± 


0.036 




.46 


+ 





.045 


17 


.07 




0, 


.02 


B 






3.55 ± 


0.036 




.49 


+ 





.045 


18 


.44 


+ 


0, 


.02 


Q0957+561 A 


0.36 


1.41 


1.49 ± 


0.085 


1 


.11 


± 





.104 


15 


.60 




0, 


.03 


B 






1.43 ± 


0.067 


\ 


.10 


+ 





.095 


15 


.68 


+ 


0, 


.03 


LBQS1009-0252 A 


0.87 


2.74 


1.75 + 


0.082 


I 


.32 


+ 





.036 


16 


.63 


+ 


0, 


.02 


B 






2.97 + 


0.089 


1 


.79 


+ 


0, 


.057 


18, 


.20 


+ 


0, 


.04 


B1030+071 A 


0.60 


1.54 


4.82 ± 


0.221 


3 


.59 


± 





.250 


15 


.88 




0, 


.02 


B 






3.83 + 


0.212 


2 


.59 


+ 





.242 


19 


.79 


+ 


0, 


.03 


HE1104-1805 A 


0.73 


2.32 


1.35 ± 


0.321 





.81 


+ 





.050 


15 


.57 


+ 


0, 


.03 


B 






1.66 ± 


0.165 





.97 


+ 





.064 


17 


.04 


+ 


0, 


.04 


PGl 1 15+080 Al 


0.31 


1.72 


1.19 ± 


0.112 





.71 


± 





.028 


15 


.71 


-j- 


0, 


.02 


A2 






1.41 ± 


0.095 





.64 


± 





.042 


16 


.21 




0, 


.03 


B 






0.69 + 


0.453 





.21 


± 





.393 


17, 


.70 


+ 


0, 


.05 


c 






1.72 + 


0.322 


1 


.14 


+ 





.342 


17, 


.23 




0, 


.04 


B 1422+231 A 


0.34 


3.62 


2.02 ± 


0.112 


1 


.47 


+ 





.082 


14, 


.41 







.02 


B 






2.16± 


0.104 


1 


.56 


± 





.076 


14 


.29 




0, 


.03 


c 






2.11 ± 


0.076 


1 


.43 


+ 





.050 


14 


.98 


+ 


0, 


.03 


D 






2.30 ± 


0.063 


1 


.54 


+ 





.045 


18, 


.14 


+ 





.02 


SBS1520+530A 


0.72 


1.86 


1.85 ± 


0.161 





.77 


+ 





.036 


17, 


.20 


+ 


0, 


.02 


B 






1.45 ± 


0.163 





.96 


+ 





.095 


18, 


.03 


+ 


0, 


.03 


B 1600+434 A 


0.41 


1.59 


2.95 ± 


0.614 


1 


.26 


± 





.170 


20 


.66 


± 


0, 


.13 


B 






1.85 ± 


0.891 





.92 


± 





.146 


20 


.47 


± 


0, 


.14 


MG2016+112A 


1.01 


3.27 


2.44 ± 


0.289 


1 


.97 


+ 





.086 


20, 


.48 


+ 





.07 


B 






3.15 + 


0.054 


2 


.12 


+ 





.071 


20, 


.50 


+ 


0, 


.05 


C 






4.23 + 


0.711 


2 


.94 


+ 





.064 


20 


.21 


+ 


0, 


.04 


B2045+265 A 


0.87 


1.28 


6.07 ± 


0.210 


3 


.14 


± 





.110 


19 


.81 


± 


0, 


.01 


B 






6.03 ± 


0.072 


3 


.22 


+ 





.081 


20 


.37 


+ 


0, 


.04 


C 






5.99 ± 


0.233 


3 


.13 


+ 





.126 


20, 


.05 


+ 


0, 


.04 


S 






5.44 ± 


0.057 


2 


.39 


± 





.072 


19 


.41 


± 


0, 


.04 


HE2149-2745 A 


0.50 


2.03 


1.30± 


0.085 





.85 


± 





.042 


15 


.67 


± 


0, 


.03 


B 






1.44 + 


0.076 





.85 


+ 





.036 


17, 


.23 


+ 





.03 


Q2237+0305A 


0.04 


1.69 


2.39 ± 


0.108 





.96 


+ 





.134 


14, 


.96 


+ 


0, 


.06 
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B 






2.01 + 


0.120 


1.75 + 


0.110 


15.46 


+ 


0.01 


C 






2.85 ± 


0.101 


1.06 + 


0.120 


15.71 


+ 


0.01 


D 






3.01 ± 


0.133 


1.39 ± 


0.050 


16.00 


± 


0.03 


QJ0158-4325 A - 


- 1 


.29 


1.63 + 


0.133 


1.34 + 


0.361 


16.47 


± 


0.03 


B 






1.64 ± 


0.173 


1.35 + 


0.331 


17.27 


+ 


0.03 


APM08279+5255 A - 


- 3, 


.87 


2.50 + 


0.036 


0.79 + 


0.054 


14.46 


+ 


0.02 


B 






4.10 ± 


0.036 


2.00 ± 


0.028 


13.64 


± 


0.02 


FBQ095 1+2635 A - 


- 1, 


.24 


1.67 ± 


0.204 


1.08 ± 


0.057 


15.62 


± 


0.04 


B 






1.33 ± 


0.076 


0.90 + 


0.042 


16.99 


+ 


0.03 


Q1017-207A - 


- 2, 


.55 


1.77 ± 


0.042 


1.26 + 


0.036 


15.66 


+ 


0.03 


B 






1.77 ± 


0.139 


1.26 + 


0.058 


17.81 


+ 


0.05 


Q1208+101A - 


- 3, 


.80 


2.51 ± 


0.124 


1.31 ± 


0.042 


15.91 


± 


0.03 


B 






2.11 ± 


0.045 


1.08 ± 


0.092 


17.55 


± 


0.02 


FBQ1633+3134A - 


- 1, 


.52 


1.77 + 


0.163 


1.07 + 


0.042 


15.78 


+ 


0.03 


B 






1.72 ± 


0.151 


1.06 + 


0.022 


17.23 


+ 


0.02 



